Low energy magnetic excitation spectrum of the unconventional ferromagnet CeRh3B2 
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The magnetic excitation spectrum of the unconventional ferromagnet CeRh3B2 was measured by 
inelastic neutron scattering on single crystal sample in the magnetically ordered and paramagnetic 
phases. The spin-wave excitation spectrum evidences high exchange interaction along the c-axis 
about two orders of magnitude higher than the ones in the basal plane of the hexagonal structure. 
Both strong out of plane and small in plane anisotropies are found. This latter point confirms that 
considering the J=5/2 multiplet alone is not adequate for describing the ground state of CeRh3B2. 
Quasielastic scattering measured above Tcurie is also strongly anisotropic between the basal plane 
and the c-axis and suggests localized magnetism. 

PACS numbers: 



INTRODUCTION 



CeRh3B2 is a Ce-based intermetallic compound with 
quite unique properties. It is ferromagnetic while most 
of the other Ce-based compounds are antiferromagnetic. 
The easy axis is in the basal plane of the hexagonal struc- 
ture and the saturation magnetization 0.4 ^b is strongly 
reduced compared to free cerium ion value (2.14 /is). 
More striking is its huge Curie temperature, Tctirie=115 
K, that is two orders of magnitude higher that what 
would be expected from applying the de Gennes scal- 
ing to GdRhgBi (Tcnrje=90 K) [l[. Different theoret- 
ical models from itinerant to localized magnetism have 
been proposed to explain the peculiarities of this com- 
ound and they are reviewed in the recent calculations 
yj. The key ingredient is certainly the very short Ce- 
Ce distance (3.08 A) along the c-axis that leads to the 
proposition of strong interatomic Ce-4/-Ce-5(i hybridiza- 
tion or even direct Ce-4/-Ce-4/ interaction ,4j. While 
initially itinerant magnetism was favoured, many exper- 
imental results point toward localized 4/ electrons with 
strong hybridization as inferred from photoemission spec- 
troscopy [5] , NMR 6] dHvA measurements [7| . Polarized 
neutron diffraction and crystal field spectroscopy give an 
evidence of a mixing of J—b/2 and J=7/2 multiplets in 
the ground state wave-function [8|, |9|. Compton scat- 
tering experiments 10|, luj show that the orbital to spin 



moment ratio is less than the expected value for Ce + ion. 
Hence the spin magnetism is enhanced in CeRh3B2 with 
respect to the orbital part and this could explain partly 
the high Curie temperature since exchange interactions 
couple spins. These overall features together with the 
strong recent interest for other / electron ferromagnetic 
systems, namely the new uranium based ferromagnetic 
superconductors [12'|, motivated us to investigate the spin 
dynamics of CeRh3B2. It is worthwhile to note that no 
trace of superconductivity was detected when ferromag- 
netism is suppressed at 7 GPa 13| in CeRh3B2 while su- 



TABLE I: Instrument configurations. The given colhmation 
is the one after the monochromator. Fixed final wave-vector 
was used with a filter on the scattered beam. PG is pyrolitic 
graphite 
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perconductivity is known to exist in Ce(Rhi_j.Ruj;)3B2 
for 0.62 < a; < 1 and ferromagnetism disappears for x > 
0.16 [3. 



EXPERIMENTAL DETAILS 

Experiments were carried out on the three axis spec- 
trometers IN12, IN22, IN20 and INS at Institut Laue 
Langevin, Grenoble. The configuration of each spectrom- 
eter is given in Table 1. For all configurations, a fixed 
final energy was used together with a focusing analyzer 
without collimation. Except for INS where measurements 
were performed around (0, 0, 2), most measurements 
were carried out around the very weak nuclear refiection 
(0, 0, 1) in order to limit the contamination by acoustic 
phonons. The single crystal was grown by the Czochral- 
ski method by using isotopically enriched ^^B (90 %) in 
order to reduce neutron absorption. The sample for neu- 
tron scattering consists of two co-aligned platelets with 
[0, 0, 1] normal. Their thickness is 2 mm in order to 
reduce the effect of Rh neutron absorption. They are of 
typical length 25 mm and width 5 mm. 
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FIG. 1: Magnetic excitation spectrum of CeRh3B2 at Q=(0, 
0, 1) for a) T = 150 K and b) 5 K. Solid lines are fit to the 
data with a quasielastic Lorentzian at 150 K and an inelastic 
Lorentzian at 5 K (See Appendix). The dashed line indicates 
the background. 



SPIN-WAVES 

Figure 1 shows magnetic excitation spectrum mea- 
sured on IN12 at Q=(l, 0, 0) for T = 5 and 150 K. 
They were fit with inelastic Lorentzian below Tcurie and 
quasielastic Lorentzian above Tcurie (See Appendix). 
These spectra give evidence of a gap in the excitation 
spectrum of about 2 meV for q=0 in the ferromagnetic 
phase (In this paper, Q=T-|-q where r is a zone center 
and Q=(Q_f/, Qif , Ql) and q=(h, k, 1)). The dispersion 
along the [1, 0, 0] direction was investigated by perform- 
ing several constant Q scans. Figure 2 shows a represen- 
tative magnetic excitation spectrum measured on IN20 
at Q=(0.1, 0, 1) for 5 and 150 K. In this work, the back- 
ground is consistently taken for each configuration as the 
flat part of the constant energy and constant Q scans. 
The dispersion along the [0, 0, 1] direction was measured 
by performing constant energy scans since the dispersion 
is much steeper in this direction than along [1, 0, 0]. 



Representative spectra are shown in Figure 3. Figure 3a 
shows measurements performed on IN20 for an energy 
transfer of 9 meV subtracted from the background and 
corrected for the Bose factor. In the ferromagnetic phase, 
the asymmetry between the two peaks located on both 
sides of (0, 0, 1) corresponds to the defocusing and fo- 
cusing conditions for the measurements. The data in the 
paramagnetic phase exhibits weak Q-dependence due to 
magnetic correlations (see next part). Figure 3b shows 
a subtraction of the data measured at 5 K and 150 K 
for an energy transfer of 40 meV on INS. For this mea- 
surement the background is not well determined due to 
optic phonon and multi-phonon contamination near (0, 
0, 2) so that the same procedure as the one shown in 
Fig. 3a is not applied. However the peak position is still 
well-defined. The resulting dispersion is shown in Figure 
4. There are three salient features : (i) the small gap at 
q=0, (ii) the linear dispersion at small q and the very 
high anisotropy of the spin- wave energy between the two 
directions. Without further analysis these facts already 
allow to draw some conclusions on the spin dynamics of 
CeRh3B2. Indeed for a ferromagnetic planar system, the 
spin-wave spectrum is linear and gapless. The observed 
gap is thus a signature of the in-plane anisotropy. There- 
fore, the spin wave is analyzed by using the following 
dispersion relation [15| : 

ul = [Ai + 25/(0) - 25/(q)] [A2 + 25/(0) - 25/(q)] 

(1) 
LOq in the energy of the mode, Ai is the axial anisotropy 
and A2 is the in-plane anisotropy. /(q) is the Fourier 
transform of the exchange integral /^ j between the mag- 
netic sites i and j with the Hamiltonian ^^ ■ — /^jS^Sj. 
Si is the spin operator at site i and S is the value of 
the spin. Table II shows the different exchange integrals 
and the corresponding distances between Ce atoms and 
the number of neighbors z. Given the shape of the dis- 
persion, the relevant exchange parameters that we have 
deduced are : Ic the nearest neighbor interaction along 
the c-axis, /i the nearest neighbor interaction in the basal 
plane and I2 , the next nearest neighbor interaction in the 
plane. Note that because measurement are performed 
only along [1, 0, 0] and [0, 0, 1], the exchange Idiag along 
[1, 0, 1] is included in Ic and h. If Idiag is strong (which 
is not really expected because the corresponding Ce-Ce 
distance is high), Ic and /i are effective parameters. We 
cannot also determine separately the exchange parame- 
ters and the energy gaps because they occurs as products 
in Eq.(l). The obtained parameters from a fit of the dis- 
persion shown in Fig. 4 are given in Table II. /i is not 
very well determined because part of I2 intervenes in the 
same way as /i in the fit. The salient feature is the large 
difference between Ic and /i (/i), the former being two 
orders of magnitude larger than the latter : Ic/Ii ~ 300 
and Ic/l2 ~ 100. The dispersion along the [0, 0, 1] di- 
rection is characteristic of the regime Ai > 2SIc and 
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FIG. 2: Magnetic excitation spectrum of CeRh3B2 at Q=(0.1, 
0, 1) for T = 5 and 150 K. Solid lines are fit to the data with 
a quasielastic Lorentzian at 150 K and a Gaussian at 5 K (See 
Appendix). The dashed line indicates the background. 



good fits are indeed only obtained for Ai > 100 meV 
that leads also to the constraint 2SIc < 20 meV. We will 
tentatively give individual estimate for Ai and Ic in the 
discussion part with extra hypotheses beyond the present 
fit. The gap at q=0 is the geometric mean of the axial 
and in-plane anisotropics: Wg=o=vAiA2 ~ 2 meV. 



PARAMAGNETIC SCATTERING 

Preliminary data concerning paramagnetic scattering 
obtained on powder sample were shown in Ref. 9]. A 
quasielastic signal was observed with a powdered average 
linewidth of 2 meV at 150 K and 5 meV at 300 K. Figure 
5 shows representative spectra taken in the present single 
crystal study at 120 K just above Tcurie for three wave- 
vectors, which establishes the Q dependence of the sig- 
nal. The data are analyzed with a quasielastic Lorentzian 
lineshape convoluted with one dimensional resolution in 
the Lo direction. This analysis gives two parameters, the 



TABLE IL Exchange integrals and anisotropies 
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FIG. 3: a) Constant energy scans performed along [0, 0, 1] 
at 9 meV at 5 K and 150 K with the background subtracted 
and corrected by the Bose factor, b) Subtraction of the raw 
data measured at 5 and 150 K by constant energy scan along 
[0, 0, 1] with energy transfer of 40 meV. Lines are fits with 
Gaussians (see Appendix). 
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FIG. 4: Spin-wave dispersion in the [1, 0, 0] and [0, 0, 1] 
directions at 5 K. Lines are fit as explained in the text. 
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FIG. 5: Energy spectra measured at 120 K. Solid lines are 
fit to the data with a quasielastic Lorentzian (See Appendix). 
The dashed line indicates the background. 



FIG. 7: The product x^fq at 120 K for q along the [1, 0, 0] 
direction. The line is a linear fit going through the origin. 
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FIG. 6: Relaxation rate in the [1, 0, 0] and [0, 0, 1] directions 
at 120 K and 150 K. Lines are fit as explained in the text ; 
The solid line corresponds to the local model and the dotted 
line corresponds to the itinerant model. The dashed line is a 
guide for the eyes. 



q-dependent susceptibility Xq ^"^^ the relaxation rate V q. 
In this work, the probed fluctuations are in-plane fluctu- 
ations since the wave- vector has predominant component 
along the c-axis (being (h, 0, 1) or (0, l-fl)) and neutron 
probes fluctuations perpendicular to Q. The obtained re- 
laxation rate is shown in Figure 6 at 120 K and for [1, 
0, 0] and [0, 0, 1] and for a few points at 150 K along 
[1, 0, 0]. The obtained values are qualitatively consistent 
with the powder averaged ones since the basal plane has 
more spectral weight than the c-axis direction. The ex- 
tent in (jj space of these excitations is smaller than that of 
the spin wave for the [0, 0, 1] direction and quite similar 
for the [1, 0, 0] direction. It seems that the relaxation 
rate extrapolates to zero for q=0. However this is not 
expected since the in-plane magnetization is not a con- 
served quantity in a planar magnet [16|. Hence a small 



finite value Vq((i = 0) may exist but our data do not 
allow us to extract it. We fit the data obtained along 
the [1, 0, 0] direction at 120 K with either T=Dh'' or 
r = ah{l + {h/ Kf') that corresponds respectively to lo- 
calized or itinerant model of ferromagnetism. Along [1, 
0, 0], we obtained D= 19(3) and z= 1.34(14) (solid line 
in Fig. 6) and respectively a=9.9(9) meV.(r.l.u.)^^ and 
K=0.6(l) r.l.u. (dotted line in Fig. 6). It is worthwhile 
to note that z is close to 3/2, the dynamical exponent 
expected for an X-Y magnet [l^|. The fit along the [0, 0, 
1] direction is not exploited due to the limited number of 
data points. Usually the ratio D/Tcurte can give some 
information on the nature of the magnetism i.e. local- 
ized versus itinerant. It is expected that < Fg > /Tcurie 
>> 1 for itinerant magnetism and that < Fg > /Tcurie 
w 1 for localized magnetism 17[ . This classification was 
applied for isotropic cubic magnets 18[ or more recently 
in orthorhombic UGe2 [I9|. It is highly questionable to 
apply it here since the value of D is quite different be- 
tween the plane and the c-axis (See Fig. 6). We could 
nevertheless notice that the powder averaged value of the 
relaxation rate [9| leads to < F^ > /Tcurie < 1 that is in 
favor of localized magnetism. This conclusion is also sus- 
tained by the value of k obtained by the itinerant model 
fit of Tq that is quite large for a temperature just above 
the Curie temperature and this analysis seems therefore 
unphysical. Hence our limited set of data obtained in the 
paramagnetic phase would suggest localized magnetism. 
Finally as compared to antiferromagnetic heavy fermion 
systems where Xg^g is constant, it is worthwhile to point 
out that here XgF^ is almost linear in q for the [1, 0, 0] 
direction at 120 K for which data were extensively taken 
(See Fig. 7). This emphases the peculiarity of ferro- 




FIG. 8: Constant energy scans performed for an energy trans- 
fer of 3 meV at 120 and 200 K along [1, 0, 0] and [0, 0, 1] 
directions. Solid lines are fit to the data with Gaussians as 
explained in the text. The dashed line indicates the back- 
ground. 



magnetic fluctuations and the almost conserved magne- 
tization. In the ferromagnetic compound UGe2, it was 
found that XgFg is linear with a large finite intercept at 
q—0 [20|. Constant energy scans performed for an en- 
ergy transfer of 3 meV are shown in Figure 8 for 120 and 
200 K for [1, ,0] and [0, 0, 1] directions. The peak 
width is very anisotropic between the basal plane where 
peaks from adjacent Brillouin zones overlap and the c- 
axis where peaks are narrow. The constant energy scans 
can be fit by Lorentzian or Gaussian lineshape and both 
fits provide almost the same width k*. This is not the 
true inverse correlation length k since the measurement is 
performed for a finite energy of 3 meV while k is obtained 
by energy integrated imaginary part of the dynamical 
spin susceptibility. The gaussian fit gives Kt g (,,=0.33(1) 
r.l.u. , K*g„^j =0.065(4) at 120 K and kJ-^qq] =0.38(5) , 

[0,0 



^,=0.100(8) at 200 K. The temperature dependence 
of this characteristic length is more important along the 
c-axis. This could explain the above mentioned fact con- 
cerning the extent in q— cj space of the fluctuations. The 
change of this extension along [0, 0, 1] between the para- 
magnetic and ordered phase could be due to a stronger 
temperature dependence of the different parameters in 
this direction. Exploiting further the paramagnetic scat- 
tering would need a detailed survey of the spin dynamics 
in the (Q,a;) space above Tcurie- 



DISCUSSION 

In this section, we give an estimate of Ic beyond the fit 
performed in section III that only gave the constraint 
2SIc < 20 meV and 2S'Ai/c « 2500. Knowing that 
Ic » Ii^h, it is tempting to describe the system as 
a set of weakly coupled ferromagnetic chains 2l|. In 
such a model Ic » Tcurie and three dimensional or- 
der occurs for Tcurie ~ 2S{S -\- 1)\/Icli for 6 neighbor- 
ing chains [22l |. However here our spin wave data show 



that 2S/c < 20 meV. Knowing that 25* ~ 1 from Comp- 
ton scattering, we deduce that Ic is not much stronger 
than Tcurie and the weakly coupled ferromagnetic chains 
model is therefore not valid here. For the purpose of an 
estimate, we therefore use the usual mean field approx- 
imation in its crude formulation, Tcurie ~ Ic, since Ic 
is the dominant coupling (In such a case, the exact for- 
mula is Tcurie = ^S{S + l)/c). This leads to the fol- 
lowing order of magnitudes : /c ~ 10 meV, Ai w 250 
meV and A2 « 0.02 meV. This is to be compared to 
the value of the out of plane anisotropy infered from the 
latest XAS results, Ai w 50 meV Q. The origin of 
the X-Y nature of the system is the strong crystal field 
anisotropy. The dispersion relation given in formula (1) 
was used for exchange dominated planar ferromagnets 
like Tb [2J| with in this case the use of total angular mo- 
mentum J instead of S*. In such a case, the anisotropy 
gaps are expressed as Ai=6Ji32 and A2=36J5i3g with 
Jn = J {J — !)...( J — n/2) and B"^ are the canonical co- 
efficients of the crystal field Hamiltonian when expressed 
with the Stevens formalism [25|. We cannot use this for- 
mulation for CeRh3B2 since for Ce^"*" with a total mo- 
mentum J=5/2 ground state, the coefficient B^ is zero. 
As stated above, the gap at q—Q is a strong indication 
that in plane anisotropy exists. For Ce^"*" this is impos- 
sible with J=5/2. Hence the multiplet J=7/2 has to 
be taken into account in the ground state wavefunction. 
This mixing between J=5/2 and J=5/2 in the ground 
state was first pointed out by the form factor measure- 
ments 8]. Concerning the anisotropy in the plane, one of 
the reported magnetization measurements at 2 K found 
a difference of 0.004 hb for the saturated magnetization 
between a and a* axis [26l |. However this result was not 
reproduced by another group that also reports possible 
sample disorientation issue in the aforementioned work 
[2 71 ] . The small value of A2 that we proposed here implies 
that its signature will only occur at very low temperature 
probably much below 2 K. 



CONCLUSION 

In the present paper, we have focused on the magnetic 
interactions in CeRh3B2 while previous neutron scatter- 
ing studies focused on single-site crystal-field contribu- 
tions. Usually rare-earth compounds are classified into 
exchange dominated or crystal field dominated systems 
[2J|. This classification is clearly not valid for CeRh3B2 
where all energy scales must be taken into account in the 
ground state and low-lying excitations: spin-orbit, crys- 
tal field and exchange. While our study does not answer 
the question of the microscopic origin of the strong Curie 
temperature, we clearly establish the huge anisotropy be- 
tween the exchange along the c-axis and in the plane as 
well as the anisotropy within the plane that precludes a 
description of the magnetism in term of J=5/2 only. The 



spin-wave spectra extending presumably up to 80 meV 
at the c-axis zone boundary has the highest energy for 
any known cerium compound. This huge value is clearly 
a combination of crystal field anisotropy (Ai) and ex- 
change (/c). The present work is aiming to stimulate 
further theoretical studies that would be able to treat on 
the same footing all the aspects (spin-orbit, crystal field 
and exchange) of the intriguing magnetism of CeRh3B2. 



APPENDIX 

In a neutron scattering experiment, the measured in- 
tensity is the convolution of the resolution function and 
the scattering function S'(Q,cj). This latter function is 
related to the imaginary part of the dynamical spin sus- 
ceptibility x" (Q, Lo) via the fluctuation-dissipation the- 
orem S{Q,uj) = (1 - ea;p(-cj/r))-ix"(Q,w). In this 
paper, we use for the paramagnetic scattering (Fig. la. 
Fig. 2, Fig. 5), a " quasielastic Lorentzian" : 



X"(q,^) 






(2) 



For the spectra measured at Q=(0, 0, 1) (Fig. lb), 
wc use the following form, that was named "inelastic 
Lorentzian" in the paper : 



x"(q,^) = ^ 



Xq^q^ 



Xq'^q^ 



(^ 



^0) 



Tl {co + ujor + ri 



(3) 



This form is phenomenologically used to reproduce the 
broadening of the peak. The origin of this anomalous 
shape is presumably the steep dispersion along the c-axis 
that is picked up by the finite resolution function. This 
effect is strongest at the minimum of the dispersion. For 
the spin wave spectra at finite q, we use a gaussian for 
5(Q,cj) both for the constant Q (Fig. 2) and constant lj 
(Fig. 3) scans, i.e for x=h,l or lu. 

S(q, uj) = S{x) = I.exp{-ln2{{x - C)/W/2f) (4) 

with / the peak intensity, C the position of the peak and 
W its Full Width at Half Maximum (in Q or cj space). 
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